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This paper deals with a topic in multivariate analysis. Consider that a sample of 
size n+1 has been collected from a p-variate normal distribution having dispersion matrix 
(o-jj'). Let aji'/n denote the usual unbiased estimate of <r,j'. Further, let 0<Cl<^u be con- 
stants such that all characteristic roots of a matrix having the Wishart distribution lie in 
the interval [I, u] with probability 1 — a. A theorem of Roy, Bose, and Gnanadesikan 
[Ann. Math. Stat. 24, 513-536 (1953); Biometrika 44, 399-410 (1957)] may be stated as 
follows: The probability is 1 — a: that every principal ininor determinant of l~^{(ijj') — {(tjj') 
and of {(Tij') — u~^{ajj') is nonnegative. The previous result may be used to prove the main 
theorem of the present paper. Theorem: The probability is at least 1 — a that the following 
system of relations hold simultaneously: u~^ajj<ajj<l~^aji; j=l, . . ., p and \<tjj'—% 
(u-^^l-^)air\<y2il-'-u-'){a^,an'y.i J9^r • 



1. Notation 

Throughout Yi={yu, . . ., Vip)' for ^ = 1, . . ., 
AT' will be a sample of size N=n-\-l from a ^-variate 
normal distribution with mean vector i={^i . . .??,)' 
and dispersion matrix 2. In short, the F, are inde- 
pendent and A^(^, 2). Except for a and X)? Greek 
letters will always refer to parameters. We will 
frequently write S=(o-y_;/) when we mean that Ojy is 
the element in the ^th row and .7'th column of 2. 
In the same spirit 2~^=(cr^^'), and A=(ajj') will be 
common notations. Here 



i = \ 



(1.1) 



2. Simultaneous Confidence Intervals 

Let 0<^l<iu be constants sucii tliat all character- 
istic roots of 8. pXp matrix having the distribution 
W(Ijn) (see [1]^ for this notation) are in the intei'val 
[/, u] with probability 1 — a. The following theorem is, 
for our purposes, a more convenient statement of 
some results due to Roy, Bose, and Gnandesikan [8,9]. 

Theorem 1. The probability is 1 — a that I'^A—i: 
and I^—u~^A are positive semi -definite. 

A proof is included since the orighnd verification 
demonstrates many other interesting results as well, 
and consequently is very indirect. There exists a 
non singular triangular matrix C such that CXC^ = I 
[1, p. 156]. Making the transformation Zi = CYi{i = 
1, . . ., N) then Z,- has the distribution N(C^, I) 

AT _ _ 

and/i^ Z) (Z,-Z)(Z,-Z)' = (^/16^Mias thedistri- 
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bution W{I,n). Tlie theorem will be proved when it 
is shown that the following three statements are 
equivalent, (i) all roots of B are in the interval 
[/, u]. (ii) B— II and uI—B are p.s.d. (positive semi- 
defuiite). (in) A—l^ and ui:—A are p.s.d. (i) and 
(ii) are clearly equivalent. That B—II and A— 11^ 
are p.s.d. together may be shown by considering 
their quadratic forms. Making the transformation 
Q=CR, where Q and R are column vectors, we 
Imve Q'{A-l:E)Q=R\a'AC-lC'i:C)R=K'(B-ir)R. 
Sinnlarly the matrices uI—B and ul^—A are p.s.d. 
together. 

An equivalent statement of theorem 1 is 

Theorem V, The probahility is 1 — a that every 
principal minor determinant of /~^^ — S and of 
ll — u~^A is nonnegative. 

Proof. The event whose probability is being calcu- 
lated in 1' is identical witli tliat of theorem 1 [3, 
theorem 46.4]. 

Clearly the previous theorem provides the means 
of determining a simultaneous confidence region for 
the elements of the dispersion matrix 2. However 
this region, call it i^, may or may not be interesting 
depending on its shape. We now begin an investi- 
gation of the^hape of .3?. Let j7'=(2:Z-iA-S is 
p.s.d.)_and .^/ -=(2:2 — i^-^ A is p.s.d.). Then i^= 
f/[\ ^ , the intersection or common part of 5^and ^. 

Theorem 2. .$^and .^are convex cones with vertices 

2==/"^ A and l,=u~^A, respectively. 

Prooj. Assume 2ie,^, 226^^ and c>0. Clearly 
l-^Ae^". Further c2i + (l-c)riAey' since I'^A- 
[c2i + (l— c)/-i^] = c(/"^A— 2i). Hence y'is a cone 
with vertex V^A. To show .i/'is convex we mus* 
demonstrate that c2i + (l—c)22€y whenever c<l 
But the quadratic form Q'\l-^A-ci:x—{^ — (^)^M- 
cQ'{l-^A-i:,)Q+{\-c)Q'{l-^A-^2)Q is nonnega 
five and 5^ is a convex cone. The proof for ;/ is 
similar. 
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According to theorem 2, .^ is the intersection of 
two convex cones. Further on in the paper (theorem 
3) we shall prove a result which implies that .^is 
bounded. Hence the confidence region ^ is convex 
and bounded. 

Now, in the case ])=2, define 

'^^1 = 0-11— '^"^^11? W2~(ri2—u~'^ai2, Ws^(r22—u~'^a22, 

di=-(l~'^—u~^)an,d2=(l~^—u~'^)ar2, 
ds={l-'-u-')a22, (2.1) 



w-~ 



\w2 wj ' 



dD- 



d, d2^ 
r/2 d^/ 



^ and ^become [Wi, W2, w^ :wi>0, Ws>0 and wl< 
Wi Wz] and [wi, 1^2, w^ : Wi <di, w^<d^ and {w2—doy < 
{wi—di){w^—d^], respectively. Solving simul- 
taneously we find that the intersection of the 
boundaries of 6^'^ and ^ takes place on the plane 
diW^-{-d2,Wi—2d2W2=didz—dl = \D\. On this plane the 
equation of the required intersection is d^wl={\D\-{- 
2d2W2 — diWs)tOs or 



dswl—2d2W2Ws+diwl— \D\ws=0 



(2.2) 



which will be an ellipse since D=(l~'*—u~^)A is 
positive definite with probability one. 

Since the vertices of ^ and ^satisfy (2.2) the ex- 
tent of W2 and u'3 in 0^ will be the same as their ex- 
tent in the ellipse. To find the extent of W2 in (2.2) 
consider W2 to be fixed and solve for w^. (2.2) is 
then a quadratic equation with discriminant 
{\D\+2d2W2y—4:didswl. If w^ is to be real, then the 
discriminant must be nonnegative ; hence the extent 
of W2 is given by the relation (\D+2d2W2y>4:didswl 
or equivalently \w2—d2/2\<y2{dids)''^. The extent 
of Ws is even easier to calculate by the same methods. 
Solving (2.2) for W2 the discriminant reduces to 
4\D\ws(ds — Ws), which yields 0<Ws<ds. Similarly 
0<Wi<di. Using theorem 1' and the eq (2.1) we 
summarize these computations as 
Lemma 1. In the case of p==2, the probability is at 
least 1—a that the following three relations hold 
simultaneously 



and 



u ^aii<(Tii<l ^aii,u ^a22<o"22<^ ^^22, 



"■12 Fi 0^12 



l- 



^ 2 



(a^a22)^ 



Now consider the general case, where p is not 
necessarily equal to 2. In theorem V we may choose 
to ignore all principal minor determinants of order 
greater than 2. If we do this, then the probability 
of the resulting event can only be increased. For 
each pair of variates, and by identical methods, we 
will obtain a system of relations just like those of 
lemma 1. 

Theorem 3. The probability is at least 1—a that 
the following S3'stem of relations hold simultaneously 

u-^ajj<ajj<l-%j,j=l,. . .,p 









3. Precision of Instruments 

A special multivariate model which frequently 
arises in connection with simultaneous measurement 
procedures [5, 11] requires that o-^^^o-^-j-o-^^ j=^, 
. . ., p and (jjy=^G^, Jt^j'- For this application, p 
becomes the number of instruments used. The 
methods of the previous section yield interesting 
results when applied to this special case. Returning 
to eq (2.2) we calculate the extent oi t^Wz—W2] this 
will yield bounds for ol in the two instrument case. 
Making the substitution Wz=t^W2, eq (2.2) becomes 

'w;2(^i+^3-2^2)+2^^2(rfi-d'2)+rfi^'-|/)|(^+t/^2)--0. 

Computations similar to those made for W2 show that 
the extent of t is \(d^—d2)±\[d^{di-{-d2,—2d2)]^. 
Equations (2.1) and theorem 1' now yield 

Lemma 2. In the two-instrument case, the proba- 
bility is at least 1 —a that the following three relations 
hold simultaneously 



cr^—an 7. 



< 2 (^11^22)-, 



2 ( J-'+u-' 
(^1— («ii— ^12) 7^ 



<^2~ (^22 — ^12) 



l-'+u~ 



2 



[an(an+a22—2av2)]% 
[(i22(an+a22—2au)V. 



For more than two instruments we may again 
choose to ignore all principal minor determinants of 
order in excess of 2. Theorem V then gives 

Theorem 4. The probability is at least 1—a that 
the following relations hold simultaneously 



max 



m 



<(7^<niin Gjj' — 
Mf L 

ax < (an-ai) - 



+u-^ , l-'-tr 



(fl'jj^. 



rr)"] 



'+M-' 



Z-' 



[fl.li («.„+«;;— 2a 



;i4 



< 



Gi<mn\ ^ [an-ai) ^ 

+ — 2 — [«'ii(«ii+^v7-2«i^.)]^ y 



plus similar inequalities involving al, 
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4. Choosing the Bounds 

The joint distribution of the characteristic roots 
of a Wishart matrix is well known, see for example 
[1, theorem 13.3.2]. Also, there is much theoretical 
work on the distributions of the extreme roots [10]. 
Tables and charts have been prepared by Pillai [7] 
and Heck [6] from which critical values of the largest 
root may be determined for certain combinations of 
n, p, and a. However, critical values of the largest 
root are of no value in this application without some 
knowledge of the smallest root. The tabulations of 
Pillai and Heck do not treat large values of their 
parameter m and consequenth^, even though there is 
a relation between the distribution of the largest and 
smallest root, their results cannot be used to deter- 
mine critical values of the smallest root in the present 
instance. The only result which seems to be ready 
for use in determining tlie bounds is an approximation 
due to P. L. Hsu (for an exposition see [2]). In our 
terminology, Hsu's result states that I (l<xl^p^2) 
is, for large n, an appi'oxiniation to the probability 
that the minimum root is at least as large as /. 

Here we present an exact treatment of the two 
variate case (^=2). Denote the roots by r and 6' 
where r<s. Our task is to determine constants 
I and u such that P(l<r<s<u) = l—a. 

P=P(l<r<s<u) 



=c 



r;, J,'., '">""«"[--?]« 



s—r)dr(ls, 



where c^^lir h^T Q) T (\^)] '• Mak 



ana' the 



transformation r = t—v, s = t+v then F/Ac may be 
written as 



t = i .7^=0 

(*u r*u-t 



V{f-lf) 2 dvdt. 



Now, integrating w^th respect to v and transforming 
the resultant expressions, we obtain (n — l)P/2c in 
the form 

, r^ {xl2Y-'e-"hlx 

/ n — y _l^ 71 — 1 W\ Hu n—l x_ 

Finally remembering the value of c and using the 
expression V7rr(r0=2^-^r C^^j T C^:r~) ^e obtain 



P{l<r<'^<u) 
=P{2l<xln<2u) 



(4.1) 



2 2 p V2y 



where xL is a chi-square variate with 2n degrees of 
freedom. 

As is usual in such problems, there is considerable 
freedom in choosing I and u so that P=l—a. In 
table l^uis taken to be + oo and in table 2, / equals 
0. The bounds of table 3 are chosen as follows: 
Determine I so that P{l<r) = l—a, then holding 
/ fixed at this value determine U2 so tliat 
P(l<r^s<U2)=^l—2a. Note that the value of 
U2 obtained in this way agrees very closely with the 
corresponding entry of table 2. It appears that the 
probabilistic dependence of the greatest and least 
root may be ignored for the purpose of this paper. 

Table 1 was used to check the accuracy of Hsu's 
approximation in the bivariate case. The approxi- 
mation runs from about 15 percent too large for 
n=()0 to about 11 ])ercent too large for n=l()0. 
Presumably the approximation would be equally 
poor in general, except for very large sample sizes. 



Table 1. Percentage points 1 of the smallest characteristic 
root of a bivariate Wishart matrix having sample size n+1 

P(Z< both roots) =1 -a 



\ 
\ 


0.005 




0.01 




0.025 




0.05 




n \ 


















\ 


















2 


1. 5815 


(-5)* 


6. 2868 


(-5) 


3. 8585 


(-4) 


1. 4999 


(-3) 


3 


5. 0125 


(-3) 


1. 0050 


(-2) 


2. 5318 


(-2) 


5.1293 


(-2) 


4 


4. 0472 


(-2) 


6. 4770 


(-2) 


1. 2161 


(-1) 


1.9810 


(-1) 


5 


1. 2641 


(-1) 


1. 8124 


(-1) 


2. 9488 


(-1) 


4. 3175 


(-1) 


6 


2. 6589 


(-1) 


3. 5734 


(-1) 


5. 3421 


(-1) 


7. 3400 


(-1) 


7 


4. 5495 


(-1) 


5. 8588 


(-1) 


8. 2809 


(-1) 


1. 0905 




8 


6.8802 


(-1) 


8. 5956 


(-1) 


1. 1672 




1. 4910 




9 


9. 5975 


(-1) 


1. 1721 




1. 5444 




1.9277 




10 


1. 2655 




1. 5183 




1. 9539 




2. 3949 




11 


1. 6012 




1, 8940 




2,3915 




2,8883 




12 


1. 9636 




2. 2957 




2. 8535 




3. 4045 




13 


2. 3498 




2. 7204 




3. 3370 




3. 9406 




14 


2. 7575 




3. 1659 




3. 8398 




4. 4944 




15 


3. 1846 




3. 6300 




4. 3599 




5. 0639 




16 


3. 6293 




4. 1108 




4. 8952 




5. 6474 




17 


4. 0906 




4. 6074 




5. 4454 




6. 2446 




18 


4. 5665 




5. 1181 




6. 0082 




6. 8533 




19 


5. 0561 




5. 6418 




6.5828 




7. 4728 




20 


5. 5586 




6. 1777 




7. 1685 




8. 1023 




22 


6. 5982 




7. 2822 




8. 3699 




9. 3884 




24 


7. 6793 




8. 4257 




9. 6068 




1. 0707 




26 


8. 7970 




9. 6038 




1. 0875 




1. 2053 




28 


9. 9472 




1. 0813 


(1) 


1. 2170 




1. 3424 




30 


1. 1126 




1. 2049 




1. 3491 




1. 4817 




35 


1.4185 




1. 5243 




1. 6885 




1.8382 




40 


1. 7376 




1. 8561 




2. 0387 




2. 2049 




45 


2. 0676 




2. 1982 




2. 3983 




2. 5797 




50 


2. 4068 




2. 5482 




2. 7655 




2. 9613 




60 


3. 1068 




3. 2698 




3. 5185 




3. 7413 




70 


3. 8313 




4. 0139 




4. 2912 




4. 5388 




80 


4. 5745 




4. 7755 




5. 0797 




5. 3504 




90 


5. 3333 




5. 5515 




5.8811 




6. 1735 




100 


6. 1051 




6. 3397 




6. 6932 




7. 0062 





*The symbol (i) means that the tabled entry is to be multiplied by 10'. 
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Table 2. Percentage points Ui of the largest characteristic root 
of a hivariate Wishart matrix having sample size n-f-l 







P(bothroots<wi) = 


1-oc 






\ a 


0.005 




0.01 




0.025 




0.05 
















2 


1.3663 


(D* 


1. 2160 


(1) 


1. 0147 


(1) 


8. 5948 


3 


1. 6162 


(1) 


1. 4568 


(1) 


1. 2416 


(1) 


1. 0740 (1) 


4 


1. 8399 


(1) 


1. 6727 


(1) 


1. 4457 


(1) 


1. 2677 (1) 


5 


2. 0476 


(1) 


1. 8735 


(1) 


1. 6360 


(1) 


1. 4489 (1) 


6 


2. 2444 


(1) 


2.0639 


(1) 


1. 8170 


(1) 


1. 6214 (1) 


7 


2. 4331 


(1) 


2. 2466 


(1) 


1. 9909 


(1) 


1. 7878 (1) 


8 


2. 6163 


(1) 


2. 4234 


(1) 


2. 1594 


(1) 


1. 9491 (1) 


9 


2. 7924 


(1) 


2. 5952 


(1) 


2.3235 


(1) 


2. 1065 (1) 


10 


2. 9654 


(1) 


2. 7631 


(1) 


2. 4840 


(1) 


2. 2606 (1) 


11 


3. 1344 


(1) 


2. 9275 


(1) 


2. 6415 


(1) 


2.4120 (1) 


12 


3. 3004 


(1) 


3. 0889 


(1) 


2. 7962 


(1) 


2.5610 (1) 


13 


3.4636 


(1) 


3. 2479 


(1) 


2. 9487 


(1) 


2. 7079 (1) 


14 


3. 6246 


(1) 


3. 4045 


(1) 


3. 0991 


(1) 


2. 8529 (1) 


15 


3. 7832 


(1) 


3. 5591 


(1) 


3. 2477 


(1) 


2.9967 (1) 


16 


3. 9400 


(1) 


3. 7119 


(1) 


3. 3948 


(1) 


3. 1387 (1) 


17 


4. 0949 


(1) 


3. 8630 


(1) 


3. 5400 


(1) 


3. 2793 (1) 


18 


4. 2482 


(1) 


4. 0126 


(1) 


3. 6842 


(1) 


3. 4188 (1) 


19 


4. 4001 


(1) 


4. 1608 


(1) 


3. 8270 


(1) 


3. 5570 (1) 


20 


4. 5505 


(1) 


4. 3078 


(1) 


3. 9694 


(1) 


3. 6944 (1) 


22 


4. 8479 


(1) 


4. 5981 


(1) 


4. 2498 


(1) 


3.9662 (1) 


24 


5. 1408 


(1) 


4. 8845 


(1) 


4. 5266 


(1) 


4. 2348 (1) 


26 


5. 4300 


(1) 


5. 1682 


(1) 


4. 8005 


(1) 


4. 5007 (1) 


28 


5. 7157 


(1) 


5. 4480 


(1) 


5. 0712 


() 


4. 7638 (1) 


30 


5. 9991 


(1) 


5. 7249 


(1) 


5. 3394 


(1) 


5. 0249 (1) 


35 


6. 6949 


(1) 


6. 4065 


(1) 


6. 0008 


(1) 


5. 6689 (1) 


40 


7. 3767 


(1) 


7. 0755 


(1) 


6. 6509 


(1) 


6.3027 (1) 


45 


8. 0476 


(1) 


7. 7342 


(1) 


7. 2922 


(1) 


6. 9284 (1) 


50 


8. 7096 


(1) 


8. 3843 


(1) 


7. 9253 


(1) 


7. 5471 (1) 


60 


1. 0009 


(2) 


9. 6623 


(1) 


9. 1689 


(1) 


8. 7586 (1) 


70 


1. 1286 


(2) 


1. 0919 


(2) 


1. 0397 


(2) 


9.9614 (1) 


80 


1. 2543 


(2) 


1. 2158 


(2) 


1. 1609 


(2) 


1. 1150 (2) 


90 


1. 3785 


(2) 


1. 3382 


(2) 


1. 2807 


(2) 


1.2327 (2) 


100 


1. 5014 


(2) 


1. 4594 


(2) 


1.3996 


(2) 


1. 3495 (2) 



*The symbol (i) means that the tabled entry is to be multiplied by 10^. 



Following a verbal presentation of the results of 
this paper, certain unpublished portions of R. 
Gnanadesikan's 1956 thesis [4] were kindly pointed 
out as being pertinent. Gnanadesikan has obtained 
general recursion formulae for calculating two-sided 
bounds on the characteristic roots of a Wishart 
matrix. In the bivariate case I have been able to 
verify, by partial integration that (4.1) is equivalent 
to Gnanadesikan ^s expression. Confidence intervals 
for the dispersion parameters were also obtained 
by Gnanadesikan. His result is different from that 
of theorem 3 and may provide an alternative 
approach to the problem. 

I express my appreciation to the National Academy 
of Sciences-National Research Council for granting 
me the opportunity to participate in their post- 
doctoral research program, and in particular to the 
National Bureau of Standards for its support of 
my research under this program. Also, I thank 
Mrs. Karen Bedeau, who wrote the code for com- 
puting tables 1, 2, and 3. 
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Table 3. Percentage points U2 to be used in conjunction with 
table 1 for obtaining simultaneous upper and lower bounds on 
the roots of a bivariate Wishart matrix 

P(Z<both roots <M2) =l-2a 



\ 










\ 2a 










\ 


0.01 




0.05 




n \ 










\ 










2 


1. 3648 


/"■IN* 


1.0066 


(1) 


3 


1. 6149 




1.2344 


(1) 


4 


1. 8387 




1. 4388 


(1) 


5 


2. 0465 




1. 6293 


(1) 


6 


2. 2433 




1.8103 


(1) 


7 


2. 4320 




1.9842 


(1) 


8 


2. 6142 




2. 1528 


(1) 


9 


2. 7913 




2.3169 


(1) 


10 


2. 9641 




2. 4773 


(1) 


11 


3. 1332 




2. 6347 


(1) 


12 


3. 2994 




2. 7894 


(1) 


13 


3. 4626 




2.9418 


(1) 


14 


3. 6234 




3. 0921 


(1) 


15 


3. 7821 




3. 2407 


(1) 


16 


3. 9389 




3. 3877 


(1) 


17 


4. 0937 




3. 5329 


(1) 


18 


4. 2470 




3. 6770 


(1) 


19 


4. 3989 




3.8198 


(1) 


20 


4. 5494 




3. 9620 


(1) 


22 


4. 8467 




4. 2424 


(1) 


24 


5. 1397 




4.5191 


(1) 


26 


5. 4288 




4. 7926 


(1) 


28 


5. 7145 




5. 0635 


(1) 


30 


5. 9978 




5.3316 


(1) 


35 


6. 6937 




5. 9927 


(1) 


40 


7. 3754 




6. 6425 


(1) 


45 


8. 0463 




7. 2831 


(1) 


50 


8. 7078 




7.9164 


(1) 


60 


1. 0009 


(2) 


9. 1633 


(I) 


70 


1. 1285 


(2) 


1.0391 


(2) 


80 


1. 2542 


(2) 


1. 1602 


(2) 


90 


1. 3784 


(2) 


1.2801 


(2) 


100 


1. 5013 


(2) 


1.3989 


(2) 



*The symbol (/) means that the tabled entry is to be multiplied by 10'. 
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